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For the one- dimensional Hubbard model subject to periodic boundary conditions we construct a 
unitary transformation between basis states so that open boundary conditions apply for the trans- 
formed Hamiltonian. Despite the fact that the one-particle and two-particle interaction matrices 
link nearest and next-nearest neighbors only, the performance of the density-matrix renormalization 
group method for the transformed Hamiltonian does not improve. Some of the new interactions act 
as independent quantum channels which generate the same level of entanglement as periodic bound- 
ary conditions in the original formulation of the Hubbard model. We provide a detailed analysis of 
these channels and show that, apart from locality of the interactions, the performance of DMRG is 
effected significantly by the number and the strength of the quantum channels which entangle the 
DMRG blocks. 
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I. INTRODUCTION 

The numerical density-matrix renormalization group 
(DMRG) method^ works best for lattice models with 
short-range interactions and open boundary conditions. 
Non-localized versions have become a major field of re- 
search, e.g., the DMRG in momentum spac o 2 i 3 i 4 and in 
quantum chemistrjaiiiiiLSaLiS. For these applications, it is 
well established that the ordering of 'lattice sites' and the 
proper choice of basis states crucially influence the con- 
vergence properties of the DMRG algorithm; for a review, 
see Refs. [Tllll2ll3ll4| . When Chan and Head-Gordon 
applied a quantum-chemistry version of the DMRG (QC- 
DMRG) to the calculation of the ground-state energy of 
selected molecules^, they found that the DMRG leads 
to significantly better results when lattice sites are re- 
ordered with the help of the Cuthill McKee algorithm 15 . 
However, using concepts inherited from quantum infor- 
mation theory, it has been shown that the Cuthill-McKee 
algorithm fails to generate an optimal ordering in gen- 
eral. In fact, it can lead to very bad configurations which 
may even prevent the DMRG algorithm from converging 
to the proper ground-state energy. 

The accuracy and convergence of the DMRG for given 
computer resources is intimately related to the entangle- 
ment of the DMRG blocks during the renormalization 
group step. Therefore, the von-Neumann entropy of the 
blocks can be used to optimize the required computa- 
tional resource s 4 ' 16 . The generation of block entropy as 
a function of system size was studied in detail by vari- 
ous groupsiii— Recently, the entropy-approach was ex- 
tended in to include the two-site entropy profile. It 
suggests a way to improve the criteria for the generation 
of basis states and a proper ordering of the corresponding 
'lattice sites'. The study of various orderings by brute- 
force algorithms confirmed the best orderings as found 
from entropy-based methods but no definite conclusions 
could be reached yel3. 

For lattice models, boundary conditions also have a 
strong influence on the performance of the DMRG algo- 



rithm. When periodic instead of open boundary condi- 
tions are used for the one-dimensional Hubbard model, 
the block entropy increases significantly with system 
size—. In order to reduce numerical efforts to solve 
problems subject to periodic boundary conditions, the 
matrix-product state description has been introduced, 
for which, however, the interaction matrices become less 
sparse, and, thus, a true gain in performance could not 
be documented yeljSi. Recent studies— 2^ indicate that 
entanglement localization and interaction localization ac- 
tually compete and should be treated on an equal foot- 
ing. The central goal remains the development of a stan- 
dard procedure to find a basis state transformation for 
a given model which minimizes the block entanglement 
and thereby optimizes the performance of the DMRG al- 
gorithm in terms of required computational resources for 
a given demand on accuracy. 

In this work we introduce a unitary basis transforma- 
tion for the one-dimensional Hubbard model with peri- 
odic boundary conditions which results in a two-chain ge- 
ometry with open boundary conditions and couplings be- 
tween nearest neighbors and next-nearest neighbors only. 
Contrary to expectation, the performance of the DMRG 
algorithm does not improve. Our analysis shows that 
the transformation opens new quantum channels which 
interfere with the kinetic-energy channel and lead to a 
substantial entanglement between the DMRG blocks. 

We organize our paper as follows. In section |H] we de- 
scribe briefly the Hubbard Hamiltonian and the unitary 
transformation to the two-chain geometry with localized 
interactions and open boundary conditions. In section lTTTl 
we discuss our numerical procedure with an emphasis on 
the control of accuracy and the data analysis. We present 
our numerical DMRG results in section llVl We find that 
the DMRG procedure is more efficient for the Hubbard 
model with periodic boundary conditions than for the 
transformed version with open boundary conditions. We 
analyze this result in terms of the influence of indepen- 
dent quantum channels in section In particular, we 
show that a super-site representation for the two-chain 
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geometry does not remedy the basic entanglement prob- 
lem of competing quantum channels. We draw our con- 
clusions in section IvTl 



B. Two-chain geometry 

Let us define the following unitary transformation for 
an even number of lattice sites, 



II. BASIS STATE TRANSFORMATION 

A. Hubbard model 

We consider the one-dimensional Hubbard model with 
uniform nearest-neighbor hopping on a finite chain of 
L s lattice sites subject to periodic boundary conditions, 

H = -tH T + UHu 

L s -1 
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where Cj a [c- a ) is the creation (annihilation) operator 

for electrons with spin a =|, J, at site j, h^ a — c^c.- CT , 
and hj = fij^ + hj^ is the occupation number at site j. 
Due to periodic boundary conditions we set cl s ,ct = co !(T . 

The single-particle interaction matrix is given by Ht- 
We use the intersite hopping parameter t as unit of en- 
ergy and set to t = 1 in the following. The two-particle 
interaction is given by Hjj and U is the strength of the 
on-site Coulomb interaction. The schematic plot of the 
model for a chain with L s = 10 lattice sites with periodic 
boundary conditions is shown in Fig. 




FIG. 1: Schematic plot of the Hubbard model with periodic 
boundary conditions for a chain with L s = 10 lattice sites. 
Solid lines denote the nearest-neighbor hopping while the on- 
site Coulomb interaction is shown by the gray shading. Num- 
bers indicate the lattice site indices. 

It is evident from Fig. ^ an d eq. Q that the single- 
particle and the two-particle interaction matrices are di- 
agonally dominated apart from the fact that Ht has two 
off-diagonal terms due to the couplings between sites 
j = and j = L s — 1 when periodic boundary condi- 
tions are employed. These two terms leads to an en- 
larged bandwidth of Ht and to a larger entanglement in 
the system as compared to the case of open boundary 
conditions. Therefore, we should find a transformation 
which reduces the bandwidth of Ht . A reordering of lat- 
tice sites cannot lead to more localized interactions as 
will be shown in section llVl Thus, we need to apply an 
appropriate unitary transformation to new basis states. 
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The back-transformation reads for j = 1,2, 
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The transformation is the result of a Lanczos basis rep- 
resentation of the kinetic energy which starts from the 
state |4>o) = c^jvac). In this way, all the operators a^ a 

are generated. The operators b^ a naturally follow as the 
antisymmetric linear combinations of the operators Cj j(T 
and c is _,- )0 .. 
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FIG. 2: Schematic plot of the transformed Hubbard model 
using the unitary transformation l p) . Single solid lines de- 
note single-particle couplings, double solid lines correspond 
to two-particle interactions. Shaded circles denote the on-site 
Coulomb interaction of strength U while empty circles cor- 
respond to strength U/2. Numbers indicate the lattice site 
indices. 

In terms of the new operators the kinetic energy be- 
comes 



H T = J2^2 



L a /2-2 
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(4) 



The geometry of the transformed model is shown sche- 
matically in Fig. [5] As seen from the figure, the trans- 
formed model displays a two-chain geometry with open 
boundary conditions. Moreover, the kinetic energy only 
couples nearest-neighbor sites of type a or b. Thus, Ht 
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is diagonally dominated as for the case of the Hubbard 
model with open boundary conditions. 

When we apply the unitary transformation to the two- 
particle interaction matrix we find 
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with the local direct, pair-hopping, and spin-flip terms 
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The schematic plot of the two-particle interaction is also 
shown in Fig. [3 It is evident that we have transformed 
the Hubbard model with periodic boundary conditions 
into a two-band problem with purely local interactions 
and nearest-neighbor electron transfers. 



C. One-dimensional representation 

The standard DMRG algorithm applies to one-band, 
i.e., single-chain geometries. By ordering the sites of 
the two-chain geometry next to each other, the Hamilto- 
nian of the transformed model takes the form as shown 
in Fig. |3 The single-particle interaction matrix con- 
tains couplings between next-nearest neighbors. The 
two-particle interaction matrices contain on-site contri- 
butions via the direct term Hu,dj and nearest-neighbor 
interactions from Hu,d> the pair-hopping term Hjj, p , and 
the spin-flip term Hjj tS . Nevertheless, all couplings re- 
main local and the model is subject to open boundary 
conditions. Therefore, we may expect that we can calcu- 
late ground-state properties more efficiently in this for- 
mulation than we can for the Hubbard model Q with 
periodic boundary conditions. 



III. NUMERICAL PROCEDURE 

A. Controlling accuracy 

As in previous wor h 8 i 16 , we use the dynamic block state 
selection (DBSS) procedure to control the numerical ac- 
curacy of the DMRG method for different models. When 
we apply the DMRG to the Hubbard model (JIJ with pe- 
riodic boundary condition and the transformed model in 
the one-dimensional geometry of Fig. [21 we fix the quan- 
tum information loss (x) which is closely related to the 
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FIG. 3: One-dimensional representation of the single-particle 
electron transfers Ht (a) and the two-particle interactions 
H\j t d (b), Hu tP (c), Hu,s (d) of the transformed Hubbard 
model. Loops in figure (b) denote on-site interactions. 



relative error of the energy of the target state. In this 
way, the block entropy as one of the most relevant DMRG 
performance parameters can be monitored for different 
model Hamiltonians. 

We choose x to make sure that the maximum num- 
ber of block states (M max ~ 3000) that our program can 
handle is not reached during the calculations, i.e., we 
do not introduce an additional quantum information loss 
besides the truncation procedure based on \- We choose 
a small minimum number of block states M m i n in order 
to make sure that its specific choice has negligible con- 
sequences and yet ensures a reliable data analysis. We 
use the entropy sum-rule as a criterion of convergence^. 
In general, five or six sweeps are carried out in order to 
make sure that the desired accuracy determined by \ has 
been reached. 



B. Performance monitoring 

A natural quantity to measure the DMRG perfor- 
mance would be the CPU time. The CPU time, how- 
ever, strongly depends on the CPU in use and a number 
of other technical issues. A software-related quantity to 
monitor is the block entropy since it determines the num- 
ber of blocks states required to reach the desired accuracy 
for the given model and, thus, the speed of the DMRG 
calculations. 

In this work we follow notations introduced in Refs. 0, 
Ih1Ih1 | . We decompose the total system into four subsys- 
tems. There are two sites, denoted by s\ and s r , with 
gi and q T degrees of freedom between the left and right 
blocks, B\ and B r , of dimensions M\ and M r , respectively. 
The blocks = B\», B-r — »B t have dimension Ml and 
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FIG. 4: Schematic plot of the system and environment block 
of DMRG. B\ and B r denote the left and right blocks of length 
I and r, and of dimension Mi and M r , respectively, where • 
stands for the intermediate sites (si and s r ) with q\ and q T 
degrees of freedom. The blocks Bl = B\» and Br = »B r 
have dimensions Ml and Mr, respectively. 



Mr, respectively. This configuration is shown in Fig. U 
The block entropies are denoted by Si, St,, S T , and 
S-r. They, as well as the site entropies S Sl and S Sl , are 
calculated form the respective reduced subsystem density 
matrices p as S — — Trplnp. The number of degrees of 
freedom per site, q\ = q T = q, is q = 4 for the Hubbard 
model and the transformed Hubbard model. 

Apart from the entropies we monitor the Schmidt num- 
ber (7) which counts the number of nonzero eigenval- 
ues of the reduced subsystem density matrix for each 
superblock partitioning, 
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where \^>t) is the wave function of the total system, 
|^ ) and |^^) arc bi-orthogonal basis states for the 
left and right blocks with the condition £\ ujf = 1. The 
Schmidt number provides information about the entan- 
glement of the subsystems when a pure target state is 
considered. In our numerical analysis we determine 7 
for a given quantum information loss \ and we demand 
LOi > 10 -15 when we determine 7. Imposing this cut- 
off value induces some minor fluctuations in the Schmidt 
number as a function of the strength of the quantum 
channels. 



C. Total quantum information 

In order to compare more rigorously the various rep- 
resentations of a quantum system, we measure the total 
quantum information, I tot , encoded in the wave function. 
To this end, we form all system blocks which contain 
Mi = 1 to Mi = L s lattice sites and sum up the quantum 
information gain of each renormalization group stepi^. If 
no truncation is applied, / t ot also equals the sum of the 
lattice-site entropies, i.e., /tot = Ylj^sj- When we use 
the DBSS approach the error in 7 to t is proportional to 
L sX - 



IV. RESULTS FROM DMRG 

Let E(L S , Nf,N±, U) denote the exact ground-state en- 
ergy of the one-dimensional Hubbard model for a finite 
chain with L s lattice sites and N a electrons with spin a 
as a function of the interaction strength U. It can be ob- 
tained from the Bethe AnsataS. In this work we study 
the paramagnetic half-filled case, N-[ = Ny = L s /2 as a 
function of U for system sizes L s < 64. All numerical 
data presented are from the results of the last DMRG 
sweep. 



A. Lattice site reordering 

A reordering of lattice sites does not effect the total 
quantum correlation in the system. Therefore, when the 
model is solved exactly, J t ot is a conserved quantity. The 
entanglement between the DMRG blocks, however, de- 
pends on the number of quantum channels in between 
them and, thus, the ordering of lattice sites has a major 
impact on the performance of the DMRG method. 
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FIG. 5: Two extreme orderings for the Hubbard chain. 

For the Hubbard model two extreme configurations are 
shown in Fig. In Fig. (a) the "communication" be- 
tween lattice sites j = and j = L s — 1 = 9 is me- 
diated by the neighboring lattice sites in between them 
(strength t) and by a direct channel (strength t'). The 
worst reordering is shown in Fig. [3] (b) where the total 
length of all communication paths between lattice sites 
j = and j = L s — 1 = 9 is maximized. In order to 
quantify these statements we consider the (total) "com- 
munication length" , 



c= E \V(j + l)-V(j)\ 



(8) 



where V permutes the L s numbers j = 0, 1, . . . , L s — 1 
into their new ordering. For example, in Fig. [S] (b) we 
set V b (0) = 0, T b (l) = L s - 1, V b {2) = 2, V b (3) = 
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L a — 2, and so on. The standard ordering in Fig. [3] (a) 
amounts to C pbc = 2(L S — 1) which can not be decreased 
by any other ordering. For the reordering in Fig. (b), 
the communication length is C mli * = 

For a more detailed analysis it is helpful to investigate 
the number of individual quantum channels between the 
left and right blocks 

[sum over all cross-links between and 
-Br, where (I + 1) is the length of £?l-] 

(9) 

We have Cf in ^(l) = 2 in the configuration of Fig. (a) 
whereas CJJJj^f (I) = L s — \2l + 2 — L s \ in the configuration 
of Fig. [5] (b). From the number of quantum channels we 
define the (total) communication length 
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which reduces to the expression (JHJ) for our example. 

More generally, not only the number by also the 
strength and the type of the individual quantum channels 
between the left and right blocks play an important role. 
In an obvious extension of (|1(J|> we may assign adjustable 
weight factors 7channoi(', U/t,...) to each channel. Typ- 
ically, the communication length C and the number of 
individual quantum channels at each link Ci n k(0 are suf- 
ficient for a first assessment of the entanglement of the 
system. The computational cost for one DMRG iteration 
step is determined by Cn n k(0 since this number does not 
depend on the ordering of the lattice sites within the two 
blocks. The overall cost of a full DMRG sweep, however, 
also depends on C due to the relationship 11U|) . In future 
applications, C may serve as a cost function to optimize 
the ordering (and the basis set). 

In Fig. we show the site entropy, the block entropy, 
and the Schmidt number from exact DMRG calculations 
for L s — 10 lattice sites at U — 1 for the two configura- 
tions of Fig. In order to make visible the differences 
in Schmidt numbers for small system sizes we include re- 
sults for x — (exact calculation) and x = 10~ 4 with 
Mmin = 4. After convergence all lattice sites posses the 
same entropy, S Si « 1.377, and /tot ~ 13.77 is the same 
for both orderings. However, the block entropy is much 
larger for the configuration [3] (b) , i.e., one needs more 
computational resources to solve the problem. Corre- 
spondingly, the Schmidt number 7(x) f° r a given quan- 
tum information loss x is larger for the worst ordering 
than for the natural ordering. In fact, the block entropy 
and the Schmidt number closely follow the number of 
cross-links between the left and right blocks, C\[ n ^(l), up 
to logarithmic corrections in the system siz o 17 i 18 . 

There is no site ordering different from the configura- 
tion (a) which reduces the number of cross-links be- 
low two, CJPj? = 2, and the communication length below 
C mm = 2(L S — 1). Thus, we conclude that the total quan- 
tum correlation in the system cannot be reduced by a re- 
ordering of the sites. Only a basis-state transformation 
might offer a way to achieve the desired entanglement 
reduction. 
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FIG. 6: Site entropy and block entropy for the half-filled Hub- 
bard model for the two the ordering criteria shown in Fig. 
for U — 1 and \ — (exact calculation). The Schmidt num- 



ber is also shown for \ 
guides to the eyes. 
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B. Basis-state transformation 

Since interactions are localized for the transformed 
Hubbard model we might expect that the entanglement 
in the system is reduced, and, thus, the problem can be 
solved more efficiently using DMRG. In Fig.[7|we plot the 
site entropy, the block entropy, and the Schmidt number 
for the non-interacting Hubbard model and the trans- 
formed model for x = 10~ 4 , M min = 64, and L B = 64 
sites as a function of the number of DMRG sweeps. Note 
that the limit U — poses a non-trivial problem for the 
position-space DMRG. 

For both models we determine the ground-state energy 
within the desired relative accuracy of better than 10~ 3 . 
The comparison of data points in Figs.0shows, however, 
that the basis-state transformation J5J did not lead to a 
significant improvement: the site and block entropies as 
well as the Schmidt number are only marginally smaller 
than for the Hubbard model with periodic boundary con- 
ditions. 

This result can again be understood from the number 
of cross links and the communication length in Fig. [3] (a). 
The number of cross links is the same for both repre- 

2, and the communication 



sentations, C^^ 



/Ttransf 
°link 



length is almost the same, too, C' transf 
C pbc = 2L S - 2, see section UTAl 
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FIG. 7: Same as Fig.E]for the Hubbard model with periodic 
boundary conditions and for the transformed Hubbard model 
with open boundary conditions for \ = 10 -4 , M m i n = 64 
(PBC) and M min = 256 (Transf.), L s = 64, and [7 = 0. 



FIG. 8: Same as Fig. for the transformed Hubbard model 
with open boundary conditions for U = 0, x = 10 -4 , and 
M m in = 64 in independent-chain geometry. 



In order to demonstrate the importance of the number 
of links and of the communication length, we treat the 
two chains for a-electrons and 5-electrons separately as 
they decouple for U = 0, see Fig. |21 For this geometry, 
the DMRG result is shown in Fig. |S| The block entropy 
profile clearly shows the absence of quantum correlations 
between the two independent chains. For this geometry 
we have C^jf chain (Z) = 1 for the a-chain and 6-chain 
separately. The communication length for both chains 
together is C two_chain = L s - 2. Therefore, it is smaller 
by a factor of two than for periodic boundary conditions, 
C*P bc = 2L S - 2. It is evident from Fig. that the maxi- 
mum of the block entropy has equally dropped by almost 
a factor of two, from 3.3 for periodic boundary conditions 
to 1.7 for the two-chain geometry, and the Schmidt num- 
ber has reduced by more than one order of magnitude. 

The situation drastically changes when the Hubbard 
interaction is switched on. As shown for U — 10 in Fig. EI 
the site entropy and the block entropy are actually larger 
for the transformed Hubbard model with open boundary 
conditions so that more block states are required to reach 
the same accuracy as in the original formulation of the 
Hubbard model with periodic boundary conditions. Note 
that the oscillation in the block entropy is related to the 
dimer configuration of the Coulomb interaction, i.e., the 
number of bonds between DMRG blocks C^ si (l) oscil- 
lates between two and eight. Therefore, the block entropy 
of the two representations is the same for every second 
RG iteration step for which there are only the two i-type 
channels between the DMRG blocks. 
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FIG. 9: Same as Fig. □ for U = 10 and x = 10' 



As seen from Fig. ^3 the total quantum information 
of the transformed Hubbard model 10) in the geome- 
try of Fig. [3 (a) is smaller than that of the Hubbard 
model with periodic boundary conditions only for very 
small values of the interaction strength, U < 0(t/L s ). 
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V. EFFECTS OF INDEPENDENT QUANTUM 
CHANNELS 

In this section we study the entanglement generation 
in more detail. To this end, we switch on perturbatively 
various coupling terms shown in Fig. [3] 
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FIG. 10: Total quantum information for the original and 
transformed Hubbard models as a function of U for various 
system sizes. For L s — 8 results are exact while for larger 
system sizes data were obtained by setting % = 10 -5 . 



Apparently the interactions of our transformed Hamil- 
tonian, albeit fairly local, generate a strong entangle- 
ment between lattice sites and subsystems because they 
and the kinetic energy act as independent and actually 
competing quantum channels. We shall investigate this 
point further in section [V] Here, we merely determine 
the communication length by adding up equally the dis- 
tances for single-particle and two-particle electron trans- 
fers. The terms Hjj^, Hjj v and Hjj iS contribute equally 
to give the estimate (7 transf w 7L s /2. The larger entan- 
glement in the transformed Hubbard model as expressed 
by ^transf > C pbc _ 2(L S - 1) implies that the DMRG al- 
locates more computational resources for the transformed 
Hamiltonian than for the Hubbard model with periodic 
boundary conditions. 

From a technical point of view, the overall CPU time 
increases by a factor of four to five also because more 
matrix multiplications are necessary. In contrast to the 
Hubbard model with periodic boundary conditions, nine 
matrix multiplications instead of two must be carried out 
during the superblock diagonalization and three times 
more operators need to be renormalized. Moreover, due 
to the new channels and the increase of entanglement in 
the system, the number of Davidson matrix multiplica- 
tion increases by a factor of two to three. 



A. Smooth interpolation between open and 
periodic boundary conditions 

First, we consider how entanglement between DMRG 
blocks is generated for the configuration shown in Fig. [S] 
where we smoothly interpolate between open boundary 
conditions (t' = 0) and periodic boundary conditions 
(t' = 1) as a function of t' for t = 1 and [7 = 0. Our 
results are shown in Fig. 1111 The total quantum informa- 
tion, the block entropy and the Schmidt number change 
smoothly as a function of t' . Thus, the second quantum 
channel opened by the periodic boundary conditions be- 
haves perturbatively. This is not always the case as seen 
in the next example. Note, however, that the perturba- 
tion leads to a small effect only as long as t' < 0(t/L s ). 
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FIG. 11: Site entropy, block entropy, and Schmidt number 
for the half-filled Hubbard model for the configuration shown 
in Fig. Q3 (a) as a function of t' for t = 1 and U — for 
i s = 16, 32 sites, and X = 10~ 4 , M min = 64. 



8 




FIG. 12: Same as Fig. [TT] for the half-filled Hubbard model 
with open boundary condition as a function of the transfer 
integral t' between next-nearest neighbors. We set t = 1 and 
U = for L E = 16, 32 and X = 10~ 4 , M min = 64. 



Electron transfer between nearest and 
next-nearest neighbors 



As seen from Fig. the site entropy, the block en- 
tropy and the Schmidt number do not change signifi- 
cantly as a function of U. Instead, they mildly decrease 
as the interaction gradually eliminates double occupan- 
cies (and holes) from the Hilbcrt space. For small in- 
teraction strengths, the block entropy and the Schmidt 
number are fairly small. This is in accord with our ob- 
servations for the Hubbard model with open boundary 
conditions. Obviously, purely local density-type interac- 
tions do not open new quantum channels and, therefore, 
they do not substantial increase the entanglement be- 
tween blocks. Density-density interactions between dif- 
ferent lattice sites behave qualitatively the same because 
they do not involve the exchange of particles between the 
blocks. Therefore, the DMRG still performs well for the 
Hubbard model with long-range density-density interac- 
tions when open boundary conditions are applied. 

The situation changes when we treat the pair-hopping 
term Hjj,p ©, shown in Fig. 0(c), together with the ki- 
netic energy H^. We ignore all other interaction terms, 
i.e., we treat H p — —H^ + UHjj.p in Fig. Again, the 
site entropy, the block entropy and the Schmidt num- 
ber decrease smoothly as a function of the interaction 
strength. In comparison with the purely local interac- 
tion Hjj,d the Schmidt number has almost doubled. 

The increase in block entropy and Schmidt number 
is very similar when we study the effect of the spin-flip 
term Hu,s in @, as shown in Fig. (d). The result of 
the analysis of H s = —Ht -f UHjj, s is shown in Fig. Hoi 
Apparently, the spin exchange between neighboring sites 



Next, we analyze the entropy generation by an addi- 
tional next-nearest neighbor electron transfer amplitude 
for the Hubbard model. The nearest-neighbor hopping is 
again set to t = 1 whereas the transfer integral between 
next-nearest neighbors t' is smoothly increased. As in the 
previous example we choose U = to avoid the effect of 
other channels. In Fig. ^| we plot the site entropy, the 
block entropy and the Schmidt number as a function of t' . 
These quantities change smoothly as a function of t' up 
to a critical value where they increase rapidly. The value 
for the rapid increase coincides with the metal-insulator 
transition point in the t-t'-U model, see, e.g., Ref. pij . 
For finite interaction strengths, the behavior of the en- 
tropies above the quantum phase transition is more com- 
plex, corresponding to the various phases of the t-t'-U 
model2i. 



C. Density-density interactions 

We now turn to the effect of the interaction terms 
in the transformed Hubbard model (0. We start with 
the analysis of the density-terms Hu,d in ©, as shown 
in Fig. |21 (b). We neglect all other interaction terms 
and keep the single-particle hopping only, as shown in 
Fig. 13(a), i.e., we analyze Hd — —Ht + UHu,d- 




2 4 6 8 10 




2 4 6 8 10 




2 4 6 8 10 

U 



FIG. 13: Same as Fig. [HI but as a function of U for the 
local-density term only, — —Ht + UHu,d, as shown in 
Fig. B|(b). 
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creates entanglement similar to the exchange of pairs. In 
comparison of the effect of Hud 011 the one-hand side and 
Hu,p, Hjj yS on the other we conclude that not only the 
number of links and their strength but also the type of 
coupling plays an important role for the entanglement. 



D. Super-site representation 

Originally, as shown in Fig. |2 the transformed Hub- 
bard model (0 is defined on a two-chain geometry with 
purely local interaction and electron transfers between 
neighboring sites. One may wonder whether the analysis 
of subsection IV CI is adequate because it is based on the 
single-chain geometry of Fig. [3J 

In order to clarify this issue, we reduce the number of 
quantum channels between the DMRG blocks by form- 
ing 'super-sites' from lattice sites of type a and b. In this 
representation a lattice with L s /2 — 1 sites and q = 16 
degrees of freedom per site is formed, plus two end sites 
with q = 4 degrees of freedom, and open boundary con- 
ditions apply. 

The two sites at the boundaries have the same site en- 
tropy for both models. In the super-site representation 
we halve the length of the system so that sites in the inte- 
rior of the chain carry an entropy which is twice as large 
as for sites in the Hubbard model with periodic bound- 
ary conditions. A comparison of the block entropies is 
only meaningful for blocks which contain the same num- 
ber of sites, i.e., I = 3, 5,7 in Fig. ^3 As expected, for 




FIG. 14: Same as Fig. [TT] but as a function of U for the 
pair-hopping term only, H p = — Ht + UHu, P , as shown in 

Fig. B| (c). 




FIG. 15: Same as Fig. [HI but as a function of U for the 
pair-hopping term only, H s — —Ht + UHu, s , as shown in 
Fig. El (d). 



-o — <r 



O PBC 
Q Supersite 



5 6 
Site index 




FIG. 16: Site entropy and block entropy for the half-filled 
Hubbard model for periodic boundary conditions shown in 
Fig. |S] (a) and for the super-site representation in Fig. |5] for 
(7=1 and \ — (exact calculation). 



these block lengths they agree for both models. The total 
quantum information is essentially the same for both con- 
figurations. This observation is again readily explained 
by the fact that there are two quantum channels between 
the DMRG blocks in both representations. This is shown 
explicitly in Fig. U3 for U = and N = 34, X = 10~ 4 
using maximum M = 350 block states. 

The advantage that the transformed Hubbard model 
in the super-site representation is only of length L s /2 + 1 
is more than compensated by the fact that during the 
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FIG. 17: Same as Fig^for U = for N = 34, X = 10" 



numerical calculation there are q = 16 degrees of freedom 
for the two intermediate sites s/ and s r in the superblock 
representation. This higher demand for computational 
resources can be reduced if instead of two intermediate 
sites in the superblock representation only a single site is 
used. The modification of the DMRG in this direction is 
possiblo 21 ! 25 . 



VI. CONCLUSIONS 

We have constructed a unitary transformation between 
basis states for the one-dimensional Hubbard model sub- 



ject to periodic boundary conditions so that open bound- 
ary conditions apply for the transformed Hamiltonian. 
Despite the fact that the one-particle and two-particle 
interaction matrices link nearest and next-nearest neigh- 
bors only, the performance of the density-matrix renor- 
malization group method for the transformed Hamilto- 
nian does not improve significantly because some of the 
new interactions act as independent quantum channels 
which generate the same level of entanglement as peri- 
odic boundary conditions in the original formulation of 
the Hubbard model. 

The total quantum correlation in the system for the 
transformed model decreases only for small interaction 
strengths and for very short chain lengths. Therefore, 
this approach cannot be used to improve the perfor- 
mance of the DMRG for reasonable system sizes. We 
have shown that the localization of interactions alone 
does not improve the performance of the DMRG. Instead, 
it is affected more significantly by the number and the 
strength of the various quantum channels between the 
DMRG blocks. 

In conclusion, our results contribute to a better un- 
derstanding of the entanglement production within the 
DMRG. We propose to implement the communication 
length for the construction of an optimal basis and the 
proper ordering of lattice sites. The expected reduction 
of the block entanglement should improve the perfor- 
mance of the DMRG for Hamiltonians with long-range 
electron transfers and interactions. 
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